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Time scales and structures of wave interaction 
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In this paper we give a general account of Wave Interaction Theory which by now consists of two 
parts: kinetic wave turbulence theory (WTT), using a statistical description of wave interactions, 
and the D-model recently introduced in Kartashova, PRE 86: 041129 (2012) describing interactions 
of distinct modes. Applying time scale analysis to weakly nonlinear wave systems modeled by the 
focusing nonlinear Schodinger equation, we give an overview of the structures appearing in Wave 
Interaction Theory, their time scales and characteristic times. We demonstrate that kinetic cascade 
and D-cascade are not competing processes but rather two processes taking place at different time 
scales, at different characteristic levels of nonlinearity and due to different physical mechanisms. 
Taking surface water waves as an example we show that energy cascades in this system occur at 
much faster characteristic times than those required by the kinetic WTT but can be described as 
D-cascades. As D-model has no special pre-requisites, it may be rewarding to re-evaluate existing 
experiments in other wave systems appearing in hydrodynamics, nonlinear optics, electrodynamics, 
plasma, convection theory, etc. 



1. Introduction. 

Wave Interaction Theory describes the interaction of 
waves in terms of energy exchange among weakly non- 
linear modes. Two main types of energy transport have 
been intensively studied: exchange of energy among a 
small number of modes, for example an isolated resonant 
triad or quartet, and time evolution of the wave field as 
a whole consisting of an infinite number of modes. 

Interaction of distinct modes. The dynamic equations 
for a resonant triad have been solved analytically by the 
end of the 19th century and for a resonant quartet 
m 2005, 0. From the 1950s on, extensive numerical and 
experimental studies of isolated triads and quartets have 
been performed, see [3[ and bibl. therein. In the 1990s 
the notion of independently evolving resonance clusters 
formed by isolated or connected resonant triads or quar- 
tets was established 0], and their kinematic properties 
investigated, [j| . Giving also a description of the dynam- 
ics of resonant clusters, the theory of nonlinear resonance 
analysis was formulated in Q. 

In 2012 the study of the interaction of distinct modes 
was extended to describe unidirectional energy trans- 
port over scales in Fourier space (this mechanism is 
called an energy cascade). The increment chain equa- 
tion method (ICEM) for computing the energy spec- 
trum of a cascade formed by distinct modes (further re- 
ferred to as D-cascade) was introduced in Q- In Q, 
the D-model of weakly nonlinear wave interaction was 
presented, which integrates in a single theoretical frame 
D-cascades and resonance clusters and gives conditions 
when which regime will take place. 

Interaction of an infinite number of modes. The first 
oceanographic studies aiming to describe wave field evo- 
lution of surface water waves in terms of energy cascades 
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TABLE I: Structures and their computation methods 



date back to the beginning of the twentieth century. 
A major breakthrough was achieved by Hasselmann in 
1962, in [9(, giving a statistical description of an energy 
cascade in terms of a wave kinetic equation. The sta- 
tionary solution of this kinetic equation (further referred 
to as K-cascade) for capillary waves has been found in 
1967 by Zakharov and Filonenko, [l(|, coining the term 
"weak turbulence". The method was extended to other 
classes of dispersive weakly nonlinear wave systems and 
renamed to wave turbulence theory (WTT), in [ilj . 

Experimental studies performed through the last 
decades revealed numerous phenomena such as depen- 
dence of the form of energy spectra on the excitation pa- 
rameters fl2j . exponential form of the energy spectrum 
[l3j . nonlinear frequencies of the modes forming a cas- 
cade [lij], etc., which could not be covered by or even 
were in contradiction to the WTT. This was addressed 
by the invention of new models of WTT, among them 
frozen WTT, mesoscopic WTT, finite-dimensional WTT, 
[HUT?]], a H °f them relying on the statistical model of a 
K-cascade. 

The paradigm changed with the advent of the D- 
cascade which uses a deterministic model for an energy 
cascade. Presently an experimentalist has the choice be- 
tween two models for an energy cascade - K-cascade and 
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D-cascade - to describe an experiment. In this Letter we 
give criteria for this choice. 

We begin with an overview of the structures regarded 
in Wave Interaction Theory - D-cascade, resonance clus- 
ters and K-cascade, their time scales and characteristic 
times. We demonstrate that a K-cascade and a D-cascade 
are not competing processes, but rather two processes, 
which take place at different time scales, at different char- 
acteristic levels of nonlinearity and which are based on 
different physical mechanisms - resonant interactions of 
linear modes (K-cascade) and modulation instability (D- 
cascade). Also the form of the cascade is different: a 
K-cascade is decaying according to power law and is in- 
dependent of the total energy in the wave system. A D- 
cascade is decaying exponentially and its form depends 
on the energy contained in the wave system. 

Taking surface water waves as an example, we show, 
that energy cascades in this system occur at much faster 
characteristic times than those required by the kinetic 
WTT but can be described as D-cascades. 

At the end of this Letter, a list of important properties 
is given which can be used to determine from a given set 
of experimental or numerical data which type of cascade 
- D-cascade or K-cascade or even both - should be used 
to describe these data. 

2. Time scales. 
Regard a nonlinear dispersive PDE of the form 

Diji(x,t) = F(ip,il> ie ,ipt,'ipxx,il>xt,'iptt,») where (1) 

Dip(x,t) = => ip(x,t) = Aexp[i(kx - 0(k)t)], (3) 
j 



D 



d_ 
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d0(k) 
dk 



^ const. (4) 



Conditions Q provide that the dispersion function 
0(fc) is a real function and not linear in k. 

For studying ([1]) with small nonlinearity we follow 
[H ] introducing a small parameter < e <C 1 and as- 
suming that the solutions of (fTJ) have again the form 
Aexp[i(kx— 0(fc)t)], but with slowly changing amplitudes 
A(te m , xe m ). Then the nonlinearity F can be rewritten 
as 

F(£ip,eip x ,eipt,eiJ) xx ,eip xt ,eiJ}tt, ••) = 



M 

E 



: m F^ (V, V*, V*, 1>** , Vte, 1>tt, ••) + 0(e M+1 ) . (5) 



The use of any multi-scale method (in essence a trans- 
formation of variables) then allows to study the proper- 
ties of a given solution to (H} at each time scale ~ l/e m 
separately. Depending on the chosen level of nonlinear- 
ity, different physical phenomena - resonance clusters, 
D-cascade and K-cascade - can occur as it will be shown 
in the Sec. 3. 

3. Structures. 
The structures described in the next subsections and are 
shown in Table U in which they are classified according 



to the physical mechanism underlying formation of each 
structure. 

3.1. D-cascade, e-lCT 1 . 
Beginning with one wave, any multi-scale method in- 
evitably yields - at time scale 1/e 2 - the nonlinear 
Schrodinger equation (NLS) in "slow" variables T, X: 



iaxA T + a 2 A X x + a 3 \A\ 2 A = 0, 



(6) 



(for simplicity written in 1+1 dimensions, for details see 
e-g- 18[)- An NLS is called focusing if 0.3/ a\ > 0, and 
therefore admits modulation instability which in turn is 
the prerequisite for a D-cascade. In the case of non- 
focusing NLS energy cascades are also observable, e.g. in 
nonlinear optics, however no theoretical model is known 
to describe them. In this section we restrict ourselves to 
the case of a focusing NLS. 

Modulation instability is the effect that, under certain 
conditions, a wave train 0q, called carrier wave, becomes 
modulated by two side-band modes with slightly differ- 
ent frequencies 0i and 02 fulfilling conditions 0i + 02 = 
20o, ki + k 2 = 2k , where 0i = O + A0, 2 = 
0o — A0, < A0 <§C 1, [19j . For each mode 0o its instabil- 
ity interval can be written and the most unstable mode 
within this interval can be computed, also for modified 
NLS-s, e.g. [13, HJ- Modulation instability is the physi- 
cal mechanism underlying the formation of a D-cascade. 

Starting with one modulationally unstable mode 0o, 
at the first D-cascade step n = 1 two side-band modes 
are generated, 0i and 0_i, and A0i is chosen in such 
a way that one of the side-bands is the most unstable 
mode within the corresponding instability interval. The 
most unstable modes are called cascading modes. A side- 
band mode, obtained at the cascade step n from the n-th 
NLS, is the carrier mode for the (n + l)-th NLS at the 
cascade step n + 1. The direct and inverse cascade are 
not symmetrical, (22l |. 

For computing the frequencies and amplitudes of the 
cascading modes, at each cascade step n the increment 
chain equation method (ICEM) is used, 0] , yielding two 
chain equations 

0n+l = 0n + n A(0 n )k n , (7) 
0-{n+l) = 0-n ~ 0- n A(0- n )k- n , (8) 

describing the formation of a unidirectional D-cascade, 
for direct and inverse cascades respectively. After n 
cascade steps, n cascading and (2™ — n) non-cascading 
modes are excited. Non-cascading modes provide spec- 
trum broadening. The chain equations may be used to 
compute the form of the energy spectrum, the direction 
of a cascade, conditions of cascade termination, etc., [8J. 
The form of D-cascade is exponential as A 2 n = p n A^, < 
p < 1 with some constant p. The form can be approx- 
imated by c\k ai + c 2 k a2 where a\,a2 are constants for 
a given dispersion function while c\ , C2 depend on the 
energy of excitation, Q- 

It follows from 0,||HJ) that the frequencies of the cas- 
cading modes are nonlinear, i.e. dependent on ampli- 
tudes. It can be shown that O,© describe exact reso- 
nances of nonlinear Stokes waves in Zakharov's equation. 
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If the nonlinearity gets too small, all modes are linear 
and a D-cascade is not formed. Instead, under some ad- 
ditional conditions [23| . resonance clusters may occur. 

Remark 1. For a D-cascade described by (0,©, the 
characteristic level of nonlinearity is e ~ 0.1 4- 0.25, and 
the chain equation is computed from the NLS. For bigger 
nonlinearity, say e ~ 0.25 4- 0.4, a modified NLS should 
be taken with nonlinear terms up to order 4; for details 
seejpol [2l| . The corresponding chain equations are given 
in p]]. 

3.2. Resonance clusters, e ~ 10~ 2 . 
In an s-wave resonance, the application of a multi-scale 
method yields a system of equations of the form (|6j with 
interconnected coefficients whose solutions satisfy reso- 
nance conditions YljLi ^(kj) = 0' Sj=i ^kf = 0. Ne- 
glecting the "slow" space variables, one gets the well- 
known dynamical systems for 3- and 4- wave resonances, 
respectively, of the form (in canonical variables Bj ) 



Wt = VlB*B 3 , iB 



V&BiB*, iB 3 = -V&BiBr, (9) 
V^A B% B3 B±, ,iB 2 — V M B^B^B^, 



^34 

{y&YBtBiB 2 



,i£ 4 = (K 3 1 4 2 )*B 3 *B 1 B 2 



(10) 



The systems (|5j). (|10p describe the primary resonance 
clusters in 3- and 4-wave systems. In these wave systems 
also independent common clusters may be formed which 
consist of primary clusters having joint modes; they form 
the bridge to distributed initial state necessary for K- 
cascade. The form of common clusters can be computed 
by g-class decomposition, 

Time-scales for 3- and 4-wave resonances are T^^ res ~ 
X/e and T4_ res ~ 1/e 2 , respectively. They are called 
dynamic time scales. 

3.3. K-cascade, e ~ 10~ 2 . 
An initial energy distribution over an infinite number of 
resonance clusters is the prerequisite of kinetic WTT, 
[llj |. I n order to go to the kinetic regime, one has to 
make assumptions (including, but not restricted to: exis- 
tence of an infinite number of primary resonance clusters 
with linear frequencies interconnected in a special way, 
a wave field close to the stationary state, homogeneity 
in space, existence of an inertial interval where energy 
is conserved) which are necessary for the mathematical 
deduction of the wave kinetic equation. 

By statistical ensemble averaging, dynamics of individ- 
ual modes (at the dynamic time scale) is averaged out, 
and field evolution is described at the much longer kinetic 
time scale. The 3-wave kinetic equation reads as 
d 



dt 



IK 



3 2 
12 



S(0 3 - 0i - 2 )<5(k 3 - ki - k 2 ) 
(BiB 2 - B1B 3 - B*B 3 )dk 1 dk 2 (11) 



and is valid at the time scale T 3 _fci„ ~ 1/e 2 . For a 4- 
wave system a similar kinetic equation is deduced at the 
time scale T^-kin ~ 1 / £ A and so on for any s-wave system 
with finite s. Time scales T 3 _fcj n , T^_f.i n , ... are called 
kinetic time scales. An s-wave K-cascade is a stationary 
solution of the s-wave kinetic equation and occurs at the 
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TABLE II: Time scales and structures of weakly nonlinear 
wave interaction 



kinetic time scale T s _km ~ l/e 2 ( s ~ 2 K K-cascade decays 
according to power law k a where a is constant for a given 
linear dispersion relation and does not depend on the 
total energy in a wave system. 

The fact that for any s-wave system the kinetic time 
scale is substantially longer than the resonance time scale 
can be understood in the following way: a wave system 
needs a long time to move from independently evolving 
finite resonance clusters via non-resonant interactions to 
a homogeneous and almost stationary state, with energy 
distributed over an infinite number of modes. 

The structures and time scales which are known in 
wave interaction are shown in Table HU To demonstrate 
the capabilities of our approach, in the next section we 
compute time scales and characteristic times for various 
water waves. 

Remark 2. Both for resonance clusters and K-cascades, 
the small parameter e must be small enough to exclude 
nonlinear frequency shift; it is usually taken as e ~ 0.01. 
If the nonlinearity gets bigger, the dispersion relation be- 
comes dependent on amplitudes yielding a nonlinear fre- 
quency shift, and neither q-class decomposition nor ki- 
netic equation method can be used. No general methods 
are known for describing resonance clusters or K-cascades 
for nonlinear frequencies. 

4. Characteristic times for water waves. 
As it is shown in Table [ill m an y s-wave system with 
s > 4 a K-cascade is formed at a much slower time scale 
than a D-cascade. The only exception is a 3-wave sys- 
tem where a D-cascade and a K-cascade are formed at 
the same time scale ~ 1/e 2 . However, as the magni- 
tudes of the small parameter e are different for these two 
cases (e r es ~ 10~ 2 and smi ~ 10 _1 ), the characteristic 
time for formation of a D-cascade is still much shorter 
than the characteristic time for formation of 3-wave K- 
cascade. We illustrate this below, computing the char- 
acteristic time for both cascades and for various types of 
water waves. 

Let us first regard surface water waves with 2 = gk 
and wave length A = Ira. Then the wave number k = 
2tt/\ « 6, 3m _1 , the wave frequency = (9, 81-6, 3) 1 / 2 ~ 
7, 86sec _1 and the wave period is t = 2n/0 m 0.8sec. 
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TABLE III: Characteristic times for formation of D-cascade, K-cascade and resonance clustering in various water waves (at 
20°C): the gravitational constant of acceleration g — 9,8m/sec 2 ; density p — 10 3 fc<?/m 3 ; the coefficient of surface tension 
a = 72, 75 ■ 10~ s kg ■ m/sec 2 . Characteristic levels of nonlinearity are taken as in the Table HT1 emi ~ MT 1 , e res ~ HP 2 



Corresponding characteristic times are: 

for a D-cascade: e ~ lO^ 1 T MI ~ t/e 2 = t/10~ 2 ~ 
80 sec, 

for 4-wave resonances: e ~ 10~ 2 =>■ T±_ res ~ t/e 2 = 
i/10~ 4 - 8000 sec « 2,5 hours, 

for a 4-wave K-cascade: £ ~ 10~ 2 T^-kin ~ £/ £4 = 
t/10 -8 ~ about 3 years. 

In Table IIIII examples are given of the characteristic 
time, this is the time, measured in wave periods, it takes 
for a wave phenomenon to occur for different types of 
water waves and typical wave lengths. Water waves have 
been chosen as an example because at least some, though 
not many, measurements of the time scale of wave field 
evolution are known. For instance, it is established exper- 
imentally that the "fast" wave field evolution of gravity 
surface waves (caused by a sudden wind increase) occurs 
after only a few dozen wave periods [24|. The difference 
between time scales, say, Tmi and T^kin is also very 
well known to physicists working on present-day wave 
forecasting systems, e.g. 1291 . 

Last not the least. In [12| the estimate is given that 
the formation of a K-cascade in surface water waves with 
wavelength of 60 meters " would require a tank of approx- 
imately 60 km" in length. As the corresponding time 
scale is of the order of years, we conclude that it may not 
be possible at all to observe this cascade experimentally: 
A tank of 60 km in length would have to be in the open 
air, and each gust of wind would destroy the assumptions 
of the kinetic WTT about the wave fields being close to 
stationarity and being homogeneous in space. 

5. Summary and conclusions. 
In this Letter we have shown that in the focusing NLS 
three structures of wave interaction may occur - D- 
cascade, resonance clusters and K-cascade. Their prop- 
erties are listed below. 

D-cascade: 
ID) nonlinearity emi ~ 10 _1 ; 
2D) generated by modulation instability; 
3D) built of nonlinear frequencies; 
4D) occurs at the dynamic time scale Tmi ~ */ £ |fji 
5D) decays exponentially; its form depends on the total 
energy in the wave system. 

K-cascade: 
IK) nonlinearity e res ~ 10~ 2 ; 
2K) generated by s— wave resonant interactions; 



3K) built of linear frequencies; 

4K) occurs at kinetic time scales, for instance, T^-kin ~ 
t/e 2 es for 3-wave interaction and T^-kin ~ V e res f° r 4- 
wave interaction; 

5K) decays according to a power law; its form does not 
depend on the total energy in the wave system. 

Resonance cluster: 
1R) nonlinearity e res ~ 10~ 2 ; 
2R) generated by s— wave resonant interactions; 
3R) built of linear frequencies; 

4R) occurs at dynamic time scales, for instance, Tz- res ~ 
t/e res for 3-wave interaction and Ti- res ~ t/e 2 es for 4- 
wave interaction. 

For any set of data, either from experiment or from 
direct numerical simulation with evolution equations, all 
of the criteria given above must be checked in order to 
decide which model - D-cascade or K-cascade - should be 
used to describe them. 

As the D-cascade is a very young concept - the incre- 
ment chain equation method has been published in March 
2012, and the D-model in October 2012, one can not ex- 
pect to find experimental studies interpreted in terms of a 
D-cascade. But, as D-model has no special pre-requisites, 
it may be rewarding to re-evaluate existing experiments, 
especially if they use high levels of nonlinearity or are 
even known to have an exponential energy spectrum, as 

in mils. 

A D-cascade as discussed in this Letter relies on the 
physical mechanism of modulation instability which is 
found "in numerous physical situations including water 
waves, plasma waves, laser beams, and electromagnetic 
transmission lines" [27|. An application area of special 
interest is the formation of extreme waves for which mod- 
ulation instability plays a central role, e.g. [28U32T ] . 
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